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Abstract. The massless Curci-Ferrari model with Nf flavours of quarks is renormalized to three 
loops in the MS scheme in an arbitrary covariant gauge with parameter a. The renormalization 
of the BRST invariant dimension two composite operator, \A^ — ac°'c^, which corresponds 
to the mass operator in the massive Curci-Ferrari model, is determined by renormalizing the 
Green’s function where the operator is inserted in a ghost two-point function. Consequently the 
anomalous dimension of the QCD Landau gauge operator, and the (gauge independent) 

photon mass anomalous dimension in QED are both deduced at three loops. 
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There has been a renewed interest recently in trying to understand the origin of confinement 
in Yang-Mills theories and QCD from the point of view of the existence of a non-zero vacuum 
expectation value of the dimension two composite operator where is the gluon field. For 
instance, lattice gauge theory studies by Boucaud et al, for example, |Q] and references therein, 
have shown evidence that this operator has a non-zero vacuum expectation value. Whilst it might 
appear that the operator itself can have no physical relevance due to its lack of gauge invariance, 
it has been argued for instance in, [^, that in the Landau gauge it is gauge invariant since it is 
one term of a more general dimension two operator which is non-local in non-Landau gauges. 
With the assumption that the non-zero vacuum expectation value (^^4“^) is present in Yang- 
Mills theories it has been the subject of various analytic and field theoretic investigations either 
using the operator product expansion, [^, ^,0, or other methods, |, I, I, I, g n, |n|, n, n. 
One interesting approach has been that of Q. There the effective potential of the composite 
operator itself is computed at two loops in SU{N) Yang-Mills theory. It is demonstrated that 
the associated auxiliary field develops a non-zero vacuum expectation value in the true vacuum. 
The classical vacuum where the vacuum expectation value remains zero is energetically unstable. 
This calculation developed the early work of [^, IT] to compute the effective potential of a 
composite operator in a field theory to two loops which is a non-trivial exercise. One of the 
reasons for this is that non-perturbative terms arise at leading order due to the presence of a 
1/(7^ term where g is the basic coupling constant. The upshot of this is that one has to compute 
several quantities to three loops including the anomalous dimension of the gluonic dimension two 
local composite operator, [^. This was achieved by the tensor correction technique developed 
in to handle massive and massless Feynman integrals to three loop order in automatic 
Feynman diagram computer programmes. For example, the method of [Q] produced a numerical 
estimate for the vacuum expectation value of in Yang-Mills theory in the context of 

this essentially perturbative probing of a non-perturbative phenomenon. Clearly those non- 
perturbative contributions from instantons are not included in this approach but this does not 
detract from its success and potential application to other situations. Indeed in this context 
there has been recent studies of similar vacuum expectation value generation problems in Yang- 
Mills theories in different gauges. For instance, in i 0, I i |T^ , [Tl| , [l^ , |l^ the condensation 
of ghost number breaking vacuum expectation values has been investigated at one loop in the 
Landau gauge as well as the more interesting maximal abelian gauge. In the latter case the 
off-diagonal ghosts gain a mass in contrast to the diagonal ghosts remaining massless. One hope 
is that the same feature occurs in the gluon sector, indicating that the centre of the group is 
special for confinement since abelian monopoles are believed to drive this mechanism. 

Given this current interest in this area it is worth noting that it can be pursued in several 
directions. Clearly a two loop extension of ||l^ , [I^, for instance, would be interesting. However, 
all the current investigations have been for Yang-Mills theories. For more realistic studies of 
the effective potential approach one needs to include Nj flavours of quarks. Therefore, the 
aim of this article is to provide the first stage in this problem which is the determination of 
the three loop anomalous dimension of the composite operator including quarks, thereby 

extending the result given in Q. To achieve this we will in fact deduce it as a consequence of 
the renormalization of a model similar to QCD in its ultraviolet properties but differing from 
it in the infrared. Known as the Curci-Ferrari model, [0] , it was believed it could be central 
to understanding massive vector bosons as an alternative to the Higgs mechanism. However, as 
it is also not a unitary model, [IK, IT, |^, it has only received renewed interest due 


to its relation to the ghost condensation problem since the Curci-Ferrari model has a feature 
similar to the maximal abelian gauge and QCD in a class of nonlinear gauges, |^, which 
is the presence of a four ghost interaction which is a crucial ingredient for the phenomenon. As 
the Curci-Ferrari model has a natural dimension two BRST invariant mass term it can be used 
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to determine the ultraviolet structure of the renormalization group functions of QCD itself since 


the mass acts as an infrared regulator, |17, 21, 22, 23, 26, 27|. The model is renormalizable 

and the renormalization group functions are in agreement with those of QCD in the Landau 
gauge. As this BRST dimension two operator is the one required for the earlier discussion we 
will determine its anomalous dimension in the Curci-Ferrari model in an arbitrary covariant 
gauge and then specialize to the Landau gauge. Whilst this may appear a roundabout method 
we can exploit several technical points to ease our computation. Moreover, we will additionally 
renormalize the Curci-Ferrari model itself at three loops and extend the lower order calculations 
of 0, m, n, m thereby addressing several problems simultaneously. 

We begin by recalling the Curci-Ferrari model which includes a BRST invariant mass term 
for the gluon. The Lagrangian is, [17|, 


L' 


CF 






c c c c 




eab fecd^J)-c„d , ajM TPkl.H _ 


( 1 ) 


where is the gluon field, c“ and c“ are the ghost and antighost fields, is the quark 
field and a is the covariant gauge fixing parameter. The covariant derivatives are given by 
-h igApT°"ip and implying = d^A^ - d^A^ - 

structure constants of the colour group whose generators are T“, 
1 < a < Na, 1 < L < and 1 < i < with Na and Np are the dimensions of the respective 
adjoint and fundamental representations. The gluon mass is m and the quark mass is expressed in 
terms of this basic scale with the parameter /3 introduced to indicate that the masses are different. 
The model is renormalizable and the renormalization group functions are known at two loops, 
1^, ^]. In the case when mP = 0 one has QCD fixed in a nonlinear gauge which unlike the Curci- 
Ferrari model is a unitary theory. The presence of the non-zero mass in (||) breaks both unitarity 
and the nilpotency of the BRST transformation, p^, p^, 20, 21, 22, 1^ 


To compute the 


renormalization constants of (jT|) we follow a two stage approach. First, we determine the basic 
renormalization constants of the fields and parameters and then deduce the renormalization of 
the mass operator treated as a composite operator inserted in a Green’s function based on the 
massless Lagrangian. This method has been used, for example, to determine the renormalization 
of the quark mass in QCD at three loops in |^, The advantage of considering the massless 
Lagrangian is that one can apply the Mincer algorithm |32|, as written in version 2.0 of the 
symbolic manipulation language Form, |^ , ^|. This is an efficient automatic Feynman diagram 
package which determines the ultraviolet structure of massless three loop two point functions 
with respect to dimensional regularization in d = 4 — 2e dimensions. More specifically to achieve 
this for each of the Green’s functions we need to examine, we generate the appropriate set of 
Feynman diagrams using the Qgraf package, [p^j, in a format which is readable in Form, [Q|. 
These are then integrated using the appropriate Mincer routine after the basic topology of the 
diagrams has been identified. As the first stage in the computation concerns the wave function 
renormalization which derive from the two point functions we note that the renormalization 
constants are defined by 


= ^A^^^ , c“ = v^c“ , cS = v^c“ , iPo = M 


go = Zgg , mo = Zmm , Oq = Z^^Zao , /?o = Z^(3 . 


( 2 ) 


where the subscript o denotes the bare quantity. We have, using the modified minimal subtrac¬ 
tion scheme. 
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where T“T“ = CpI, Tr = Tph“^, ^“cdjfecd _ -g Riemann zeta function 

and a = g^/(167r^). We have checked our routines and programmes which determine these 
expressions by first running the files for the usual QCD Lagrangian as input before replacing the 
file defining the Feynman rules with those for the Curci-Ferrari model. We found exact agreement 
with all known QCD expressions for arbitrary a, including the more recently determined ghost 
anomalous dimension, |36, ^]. In all our calculations of the renormalization constants 

we follow the technique of |^| by computing the Green’s function in terms of bare parameters 
and then rescaling them at the end via @ which is equivalent to the subtraction procedure. 
The remaining infinities are removed by fixing the associated renormalization constant for that 
particular Green’s function. 


To proceed further we need to confirm that with these wave function anomalous dimensions, 
the correct coupling constant renormalization emerges in the Gurci-Ferrari model. As this ought 
to be gauge independent for all a this will provide a stringent check on (^) to @. However, 
to achieve this one must consider several of the three point vertices of (|l]) which, to apply 
the Mincer algorithm, requires an external leg momentum to be nullified. In this case there 
is the potential difficulty that a spurious infrared infinity can emerge in the answer due to 
infrared infinities at the nullified vertex. Therefore, it is appropriate to either consider those 
Green’s functions where this problem never arises in the first place or introduce the infrared 
rearrangement procedure, ^], which is difficult to automate. We have chosen the former 
which can be made automatic for the coupling constant renormalization but note that the 
method will also be central in deducing the gluon mass renormalization. We have computed 
the quark-gluon and ghost-gluon vertex renormalization at three loops and found that for both 
QGD and the Curci-Ferrari model the same gauge invariant coupling constant renormalization 
constant emerged which agrees with the original result of 1^^. We found 
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This together with (P) means the basic renormalization group functions for the massless Curci- 
Ferrari model at three loops are 


7a(o) = 


[(3a - 13 )Ca + STpNf] - 


+ 


+ 


^a^ + 11a — 59^ C\ + AOCaTpNj- + 32CpTpNj- 

(bia^ + 909a2 + (6012 + 864C(3))a + 648C(3) - 3986o) 

- (2304a + 20736C(3) - 58304) C\TpNf + (27648C(3) + 320) CACpTpNf 


- 9728CaT^N^ - 2304ClTpNf - 5632CpT^N^ 


1152 


+ Oia^) 


la{a) = - [(3a - 26)17^ + lOTpiV/] 


12 


a 

J 16 


^a^ + 17a — 118^ C\ + 8QCaTpNj + GACpTpNf 

(27a^ + 558a2 + (4203 + 864C(3))a + (648C(3) - 39860)) 

- (1224a + 20736C(3) - 58304) C\TpNf + (27648C(3) + 320) CACpTpNf 


- 9728CATpNf - 239AClTpNf - 5Q32CpT^N^ 


1152 


+ O(a^) 


7c (a) = 


(a — 3)Ca— + ^3a^ — 3a — 95j Ca + AQCATpNf 


0 


48 


+ 


(l62a^ + 14850^ + (3672 - 2592C(3))a - (1944C(3) + 63268)) C\ 

- (6048a - 62208C(3) - 6208) C\TpNf - (82944C(3) - 77760) CACpTpNf 


f5{a) = - 


+ 9216CaT^N^ 

fc, - iT.iV, 


a 


J 6912 


- 


+ O(a^) 

y Cl - 4CpTpNf - ^CATpNf 


+ 


2830CiTpNf - 2857C1 + 1230CACFrFlV/ - 316CAr|lV/ 

-I 

- lOSC^TpNf - 264CpT^N/j — + O(a^) . 


( 8 ) 


It is worth noting that as a consequence we have verified the three loop gauge independence of 
the /3-function by computing in this particular nonlinear covariant gauge. We have chosen not 
to re-compute the quark mass anomalous dimension since from the two loop calculation, [ ^ , 
it is clear that the same three loop gauge independent expression as ||3^, will emerge. For 
completeness, we note the anomalous dimension of the gauge parameter, a, in our conventions, 
(1^, is 


7A(a) + 7a (a) = a 
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which implies that the Landau gauge remains as a fixed point of this renormalization group 
function at three loops. 

Armed with these basic renormalization constants we have deduced the mass renormalization 
by considering the composite operator 
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( 10 ) 


in the massless Curci-Ferrari model and renormalizing it by inserting it in an appropriate two 
point function. In it was verified that this operator is multiplicatively renormalizable at two 
loops where the one loop check was established in |^. We expect this is an all orders property 
intimately related to the fact that it is BRST invariant. Indeed the interplay of renormalizability 
and BRST invariance of this operator has been explored at two loops in [44|. Clearly, we need 
to be careful which two point function O is inserted into, due to the problems noted earlier. 
There is an additional potential difficulty in this case in that the operator must not be inserted 
at zero momentum. In other words a momentum must flow through the operator, otherwise an 
incorrect result could be obtained for the operator anomalous dimension. An excellent example 
of such pitfalls has been elegantly expounded in the context of the one loop renormalization of 
in QCD in |45]. The upshot is that with O inserted in a two point function with a non¬ 


zero momentum that Green’s function is in fact a three point function. To apply the Mincer 
algorithm an external momentum must be nullified which clearly cannot be that passing through 
the operator. Instead it must be the external momentum associated with the field in which the 
operator is inserted. To determine the three loop result the only possibilities are the gluon and 
ghost fields. Inserting in a quark two point function would require a four loop calculation due 
to the absence of a tree term. Gluon external legs would inevitably lead to an infrared problem 
even at one loop which we are trying to a avoid so we are forced us to consider a ghost two 
point function. It transpires, by considering the way the momentum is nullified in this case, that 
spurious infrared infinities cannot arise, in much the same way that they do not in the three 
point function evaluations for the earlier coupling constant renormalizations. Hence, we have 
renormalized (c“(pi)0(p3)c^(p2)) with Mincer using this procedure where Pi + P 2 + Ps = 0 
and Pi = 0. Allowing for the ghost wave function renormalization of the external fields we find 
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Clearly the one and two loop terms agree with the known results for the mass renormalization 
in the Curci-Ferrari model itself, |Q, 21, 28, Moreover, the three loop answer is derived 

using the same converter routines used for the coupling constant renormalization. Also it is 
straightforward to check that like (|^ to (^) the triple and double pole terms with respect to e 
at three loops can be predicted from the one and two loop terms and these values in ( |TT|) are in 
exact agreement for all a. Hence we find 
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The only other check on this result is the Yang-Mills expression computed in for SU{Nc} 
and a = 0 with the tensor correction method. Unfortunately we do not find agreement with 
1^] for two reasons. First, the term of (|l^) involving 302428 is 377452 from equation (24) of |||] 
when compared over the same denominator. (The term 1944^(3) is in agreement with that given 
in Q.) However, it is implicit in ^ and explicit in [|^] that to deduce the corresponding Zo 
an incorrect value of the three loop Landau gauge gluon anomalous dimension has been used. 
Allowing for the correct value, (^), the expressions still differ by the quantity 19/24. We are 
confident, however, that our result (^) is in fact correct for the following simple reason, which 
to our knowledge has not been noted before. If one considers for the moment the one and two 
loop Yang-Mills results for a = 0 it transpires that the coefficients are given by the sum of the 
gluon and ghost anomalous dimensions in the Landau gauge. Indeed it is also apparent that this 
holds for Nf ^ 0. Turning to the three loop expression, (13), it is straightforward to observe 
that the same feature emerges there too. Of course such a property could be accidental but it 


justifies our confidence in the veracity of (13). Moreover, it would also suggest the existence of 
some underlying Slavnov-Taylor identity. If so then Zq is not an independent renormalization 
constant and this would reduce the number of such constants required to render (ED finite. 
The search for and construction of such a Slavnov-Taylor identity is beyond the scope of this 
article but clearly its prior existence would have rendered a certain amount of our computation 
unnecessary. However, if constructed it might play a crucial role in simplifying the calculation 
of similar anomalous dimensions in other gauges such as the maximal abelian gauge. 

Having deduced (El) in the Curci-Ferrari model we can determine the Landau gauge value 
which will correspond to the QCD value in the same gauge. We found 
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which also corresponds to the mass renormalization in the Curci-Ferrari model, , at three 
loops and together with (|^ completes the full three loop renormalization. One consequence of 
our calculation is that we can quote the value for the operator O in QED. Setting Tp = Cp = 1 
and Ca = 0 we find, from (p^). 


70 (o) I QED 
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for all a. 


To conclude with we note that the new expression for ( p!3[ ) will alter the two loop predictions 
made in Q for the numerical estimate of the operator vacuum expectation value generated by 
the effective potential method. Whilst not detracting from the achievement of that tour-de-force 
it would be interesting to repeat those calculations to explore the effect the inclusion of quarks 
has on the vacuum expectation value estimates for the more realistic case of QCD. This will 
require (^) but also needs the full two loop effective potential of the composite operator with 
quarks computed in the Landau gauge. 
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